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8: 

£Nj ■ A superconductor is influenced by an applied magnetic field. Close to the transition temperature 

T c fluctuations dominate and the correlation length £ increases strongly when T c is approached. 
However, for nonzero magnetic field H there is another length scale Lh = \l $o / a H where a is a 
universal amplitude. It is comparable to the average distance between vortex lines. We show that 
the correlation length is bounded by this length scale, so that £ cannot grow beyond Lh- This 
implies that type II superconductors in a magnetic field do not undergo a phase transition to a 
state with zero resistance. We sketch the scaling theory of the resulting magnetic field induced finite 
size effect. In contrast to its inhomogeneity induced counterpart, the magnetic length scale can be 
varied continuously in terms of the magnetic field strength. This opens the possibility to assess the 
importance of fluctuations, to extract critical point properties of the homogeneous system and to 
derive a lower bound for the length scale of inhomogeneities which affect thermodynamic properties. 
Our analysis of specific heat data for under- and optimally doped YBa2Cu307_5, MgB2, 2H-NbSe2 
and Nb77Zr23 confirms this expectation. The resulting lower bounds are in the range from 182 A 
to 818 A. Since the available data does not extend to low fields, much larger values are conceivable. 
This raises serious doubts on the relevance of the nanoscale spatial variations in the electronic 
characteristics observed with scanning tunnelling microscopy. 
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I. INTRODUCTION 



A primary effect of a magnetic field on a type II superconductor is to shift the specific heat peak to a lower 
temperature. Indeed, this shift provided some of the earliest evidence for type II superconductivity [1]. In a mean- 
■ field treatment this behavior can be comprehended as follows: To retain the continuous character of the zero field 
transition the Ginzburg-Landau correlation length £ = £o |1 — T/T c \ is supposed to diverge along the so called 
t*^ upper critical line of continuous phase transitions defined by 

^ : tc2{ H)=i-^m = ^ tf c2 (T)=Afi--V a) 



Calculations of the specific heat within the framework of the Abrikosov theory and its generalizations [2] yield 
at this mixed state to normal state phase boundary for the jump in the specific heat coefficient the expression 



,2 



, where the factor of proportionality is of order unity and takes the effect of the 

T c2 



Ac{H c2 ,T c2 )/T c2 <x (dH c2 /dT) 
(— i ■ vortex lattice into account. Matching with the zero field jump then requires that 

8 ' Ac(H c2 ,T c2 )/T c2 T c2 (H) = _ 2^ 2 ff 

^ 1 Ac(T c )/T c T c $q 1 ' 



Thus, the magnetic field reduces both, the transition temperature and the jump in the specific heat. Although 
this mean-field scenario describes the experimental data of various type II superconductors with comparatively large 
correlation length amplitude £ , including MgB 2 [3], 2H-NbSe 2 [4], and Nb7yZr 23 [5], reasonably well, it neglects an 
essential effect of an applied magnetic field. In analogy to liquid 4 He in a uniformly rotating container [6], the magnetic 
field applied to a type II superconductor is an external perturbation which drives the system away from criticality. 
Close to the transition temperature T c the correlation length £ increases strongly when T c is approached. However, 
for nonzero rotation frequency f2 or magnetic field H the average distance between vortex lines , Ln rs y (ttH) /uia^I 
or Lh ~ \/&o/ (H) is a limiting length scale. For this reason £ cannot grow beyond Lq^h- Consequently, for finite £1 
or H, the thermodynamic quantities like the specific heat are smooth functions of temperature. As a remnant of the 
singularities along T c2 (H) the quantities exhibit a maximum or an inflection point along T p (H) . Such a finite size 
effect also occurs in zero magnetic field due to inhomogeneities with length scale L. 
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In this context it is instructive to refer to the early literature focused on fluctuation effects in type II superconduc- 
tors. Implicit in most of these treatments is the assumption that fluctuations do not interact; that is, only Gaussian 
fluctuations are considered [7]. Although this assumption breaks down in the critical region, it captures some of the 
qualitative aspects of fluctuations in zero magnetic field. If the Gaussian approximation is used to calculate ther- 
modynamic properties of a type II superconductor near the mean- field phase boundary T C 2 (£T), singular behavior of 
thermodynamic and transport properties is predicted to occur [7]. However, the fluctuations of a bulk superconductor 
in a magnetic field become effectively one dimensional, as noted by Lee and Shenoy [7]; and a bulk superconductor 
behaves like an array of rods parallel to the magnetic field with diameter Lh ■ Since fluctuations become more impor- 
tant with reduced dimensionality and there is a limiting magnetic length scale Lh it becomes clear that the critical 
line T C 2 (H) is an artefact of the approximations. Indeed, calculations of the specific heat in a magnetic field which 
treat the interaction terms in the Hartree approximation and extensions thereof, find that the specific heat is smooth 
through the mean-field transition temperature T c2 (H) [8-10]. In the context of finite size scaling this is simply due 
to the fact that the correlation length of fluctuations which are transverse to the applied magnetic field are bounded 
by the magnetic length Lh cx H^ 1 / 2 . 

A main point of this paper is to show that in a type II superconductor the correlation length cannot grow beyond 
the magnetic length scale 




where a is a universal amplitude when critical fluctuations dominate. Hence, the correlation length £~ (t) is bounded 

by 

C(t)=&\t\- V = L H , (4) 

where v is the critical exponent of the correlation length. Beyond the mean-field approximation it differs from v = 1/2. 
Due to this finite size effect the maximum of the specific heat peak is shifted to a lower temperature T p given by 

Although this finite size scaling result agrees formally with the mean-field expression for t C 2 (Eq.(l)) with v = 1/2 
and a = 2n, t p (H) is not a line of continuous phase transitions. Along this line the specific heat peak adopts its 
maximum value and the correlation length attains the limiting length scale Lh ■ On the other hand, there are sample 
inhomogeneities with length scale L, as well. They imply that the correlation length £~ cannot grow beyond L, 

r (t) = l. (6) 

Thus, the zero- field specific heat singularity is rounded. Replacing in Eq.(5) Lh by L, it is seen that the maximum of 
the zero field specific heat shifts to the lower temperature T p (L). However, in contrast to the magnetic counterpart, 
where the length scale Lh oc H~ x I 2 can be varied at will in terms of the magnetic field strength, the length scale of 
the inhomogeneities remains fixed. In type II superconductors where inhomogeneities are unavoidable, two limiting 
regimes, characterized by 

L H <L, L h >L, (7) 

can be distinguished. For Lh < L the magnetic field induced finite size effect limits £ to grow beyond Ljj, while for 
Lh > L it is the length scale of the inhomogeneities. Since Lh can be tuned by the strength of the applied magnetic 
field (Eq.(3)), both limits are experimentally accessible. Lh < L is satisfied for sufficiently high and Lh > L for low 
magnetic fields. Thus, the occurrence of a magnetic field induced finite size effect requires that the magnetic field and 
the length scale of inhomogeneities affecting thermodynamic properties satisfy the lower bounds 

It is important to recognize that these bounds hold, whenever fluctuations are at work. Noting that thermodynamic 
and transport measurements are usually performed up to 12 T, it becomes clear that the magnetic field induced finite 
size effect allows to trace fluctuations, to extract critical properties of the homogeneous system and to derive a lower 
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bound for the length scale of inhomogeneities which affect the thermodynamic properties of type II superconductors. 
Such an analysis may also help to distinguish between possible mechanisms where the length scale of the inhomo- 
geneities enters in an essential way. Below, we give a summary of the main results of this paper. A related analysis 
of expansivity measurements on YBa2Cu307_,5 exposed to a magnetic field was recently performed by Lortz [11], in 
terms of a comparison with a model for uniformly rotating He [6] . 

We show that in isotropic extreme type II superconductors and nonzero magnetic field H the correlation length 
cannot grow beyond the limiting length scale Lh = yf&a] JaH), where a is a universal amplitude. Since a turns 
out to be of order one, Lh is comparable to the average distance between vortex lines. Given this limiting length 
scale, we derive the scaling expression for the line T p (H) , along which the correlation length attains the limiting 
magnetic length scale Lh- It turns out that close to zero field criticality the ratio between T p (H) and the melting line 
T m (H) is a universal number. These results, extended to anisotropic type II superconductors are found to be in good 
agreement with the experimental results of Schilling et al. [12] and Roulin et al. [13] for YBa2Cu307_,5. Moreover, 
we find that the data of Roulin et al. [13] for the temperature and magnetic field dependent specific heat coefficient, 
properly rescaled, falls on a single universal curve, which is the finite size scaling function. The phase diagram which 
emerges consists of the lines T p (Hi) and T m (Hi), where H t is the i-th component of the applied magnetic field H. 
Our analysis of the data for the optimally doped sample [13] provides consistent evidence for critical properties falling 
into 3D-XY universality class [14-16], while the data of the underdoped sample [17] uncovers the characteristics 
of a 3D-2D crossover [18,19]. Although the mean-field scenario describes the experimental data of various type 
II superconductors with comparatively large correlation length amplitudes, including MgB2 [3], 2H-NbSc2 [4] and 
Nb77Zr23 [5], remarkably well, we find that this treatment underestimates the relevance of fluctuations. Indeed the 
experimental data for the specific heat coefficient is fully consistent with a magnetic field induced finite size effect, 
although the critical regime is not attained. As a by-product we deduce for the length scale of the inhomogeneities 
affecting the thermodynamic properties of YBa 2 Cu 3 07_5, MgB 2 , 2H-NbSe2 and Nb77Zr 2 3 lower bounds ranging 
from 182 A to 814 A. Surprisingly enough, 814 A applies to the cubic superconducting alloy Nb77Zr 2 3- Since the 
available data does not extend to very low fields larger values are conceivable. In any case, these lower bounds 
are consistent with the estimates obtained from a finite size scaling analysis of the zero field specific heat data of 
nearly optimally doped YBa2Cu307_,5 single crystals, giving length scales ranging from 290 to 419 A [18,20]. This 
raises serious doubts on the relevance of the nanoscale spatial variations in the electronic characteristics observed in 
underdoped Bi-2212 with scanning tunnelling microscopy [21-24]. As STM is a surface probe, it appears to be that 
these nanoscale inhomogeneities represent a surface property. Furthermore, in spite of the consensus that in type II 
superconductors in an applied magnetic field a phase transition to a state of zero resistance should occur, we show 
that the magnetic finite size effect points to the opposite conclusion. This is in agreement with the work reported in 
Refs. [25] and [26]. 

The paper is organized as follows: In Sec. II we present a short sketch of the finite size scaling theory [27,28] and show 
that in isotropic extreme type II superconductors and nonzero magnetic field H the correlation length cannot grow 
beyond the limiting length scale Lh = \/$o/ ( a H), where a is a universal constant. Since a turns out to be of order one, 
Lh is comparable to the average distance between vortex lines. Given this magnetic field induced limiting length scale, 
we derive the scaling expression for the line T p (H), along which the correlation length attains the limiting magnetic 
length scale Lh ■ It turns out that close to zero field criticality the ratio between T p (H) and the melting line T m (H) 
is a universal number. These results arc then extended to anisotropic type II superconductors. Moreover, we find 
that the data for temperature and magnetic field dependent specific heat coefficient, properly rescaled, should fall on 
a single universal curve, which corresponds to the finite size scaling function. In Sec. Ill we analyze experimental data 
for the temperature dependence of the specific heat coefficient in various magnetic fields in terms of the inhomogeneity 
and magnetic field induced finite size effects. We concentrate on under- and optimally doped YBa2Cu307_,5, where 
single crystal data for fields applied parallel and perpendicular to the ab-plane are available [12,13,17] and extend to 
the critical regime. Accordingly, the full potential of the magnetic field induced finite size effect can be explored. The 
analysis of the data for the optimally doped sample provides consistent evidence for critical properties falling into 
3D-XY universality class, while the data of the underdoped sample uncovers the characteristics of a 3D-2D crossover 
[18,19]. Although the mean-field scenario describes the experimental data of various type II superconductors with 
comparatively large correlation length amplitudes, including MgB 2 [3], 2H-NbSe2 [4] and Nb77Zr 2 3 [5], remarkably 
well, we find that this treatment underestimates the relevance of fluctuations. Indeed, the experimental data for the 
specific heat coefficient is fully consistent with a magnetic field induced finite size effect, although the critical regime 
is not attained. As a by-product we obtain for the length scale of the inhomogeneities affecting the thermodynamic 
properties of YBa2Cu 3 07_,5, MgB 2 , 2H-NbSe2 and Nb77Zr 2 3 lower bounds ranging from 182 A to 814 A. Surprisingly 
enough, 814 A applies to the cubic superconducting alloy Nb77Zr 2 3. Since the available data does not extend to very 
low fields larger values are conceivable. This raises serious doubts on the relevance of the nanoscale spatial variations in 
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the electronic characteristics observed in underdoped Bi-2212 with scanning tunnelling microscopy [21 24]. Although 
we concentrate on the specific heat, it is shown that the magnetic field induced finite size effect also affects the other 
thermodynamic and the transport properties, including the magnetoconductivity. 



II. SKETCH OF THE SCALING THEORY FOR THE MAGNETIC FIELD INDUCED FINITE SIZE 

EFFECT 

The finite size scaling theory is based on the assumption that the system feels its finite size when the correlation 
length becomes of the order of the confining length L. For a physical quantity Q this statement can be expressed as 
[27,28] 

Q(t,L) = Q(t,L = oo) f(x), (9) 

where 

x = Trrr — v t = ^F~ l - ( 10 ) 

£(t,L = oo) T c 

L denotes the relevant confining length, T c the transition temperature and £(t,L = oo) the correlation length of 
the bulk system. Thus, the correlation length, (t, L = oo) =£=•== £^ \t\~ v , with critical amplitude critical 
exponent v, where ± = sign(t), cannot grow beyond L as t — > 0. Thus at T = T p , where 

^ (T P ) = L, (11) 

there is no singularity and the transition is rounded. Due to this finite size effect the maximum of the specific heat 
peak occurs at a temperature T p shifted by 



\t P \ 



so 



(12) 



Given the singular behavior of the specific heat coefficient of the perfect system 

% = — |*r a + B ± = UU V 
1 a 

the magnitude of the peak located at T p scales then as 



r = — \ t r a + B± = -[7±) \r a +s ± , as) 



^> = ^ M -^ +B - = Lj^i(^y\ B ,. (14) 

l p a a V c \ L J 

Here we used the universal relation [14,15] 

A ± V c =(if t ) s > V? = {£f (15) 

between the critical amplitudes of the specific heat coefficient A ± , correlation length and correlation volume V c . a is 
the critical exponent of the specific heat singularity and are universal numbers. Since superconductors fall in the 
experimentally accessible critical regime into the 3D-XY universality class, with known critical exponents and critical 
amplitude combinations, we take the universal properties, including [16] 

A + 

a = 2-'3v = -0.013, v = 0.671, — = 1.07, R~ = 0.815, (16) 

for granted. In superconductors the order parameter is a complex scalar for which one has to define below T c the 
transverse correlation length £ T . It is proportional to the hclicity modulus, related to the London penetration depth 
via [14,15,18] 



167r 3 fc B rA 



. (17) 

'o 
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Since real systems are not only inhomogcneous but are almost always impure, it is important to clarify whether 
or not quenched disorder affects the critical behavior. In general, the critical behavior of disordered systems is 
complicated, but there is a criterion, due to Harris [31]. It states that the critical behavior of disordered systems does 
not differ from that of the pure system when a < 0. This is the case when the critical exponent a of the specific heat 
is negative. 

Rewriting Eq.(14) in the form 

we realize that the L dependence of the specific heat coefficient at T p allows to determine critical amplitude and 
exponent combinations. A useful generalization of this relation is obtained from Eqs.(9) and (10), namely [32,33] 

^-^/(£Hr)-. W ..-^-(ff-.(|f. 

where g (y) is the finite size scaling function of its argument. In the limit L — > oo and at small but fixed values of t, 
g (±y = oo) = 1. When we leave L fixed and consider the limit ±t — > 0, we obtain g (y) oc \y\ a , while at T p , where 

V = Up = tp {L/Cq) 1 ^ = — 1, Eq.(14) requires that g {— 1) = 1. This scaling form allows to determine how far the 
specific heat of a system with restricted extent can be described in terms of a universal scaling function. 

The study of finite size effects in the A-transition of 4 He, as the bulk is confined more and more tightly in one 
or more dimensions [32,33], as well as computer simulations of model systems [28], have been a very good ground 
for testing such finite size scaling predictions. In this context it is essential to recognize that a finite size scaling 
analysis allows to assess the importance of critical fluctuations and to extract universal critical point properties of the 
homogeneous system. Moreover, a systematic analysis also allows to determine the length scale and even the shape 
of the inhomogeneities [32,33]. 

In type II superconductors and nonzero magnetic field H there is an additional limiting length scale. To derive this 
magnetic length scale we note that the free energy density of an anisotropic type II superconductor adopts below but 
close to the zero field transition temperature T c the scaling form [18,30,34,35] 

f s = j^G^(z), G 3 -(0) = l,z=± Ue x H x f + (e y H y f + (e z H z ) 2 ) 112 . (20) 

?y ?z 

is the transverse correlation length along direction i, diverging in zero field as = £,J_o\t\~ v , Q% is a universal 
number [14,15] and G$(z) a universal scaling function of its argument. Furthermore, ± = sign(t), t = T/T c — 1 
and H = (H x ,H y ,H z ) is the applied magnetic field. In zero field (i.e. z = 0) one recovers for G~(0) = 1 the 
scaling expression of the 3D XY universality class [14,15], leading to the universal relation (15). In the isotropic case 
(£j = £j = £j) this scaling form was confirmed by perturbation theory [36]. This scaling should give an adequate 
description of the physics of extreme type II single crystal superconductors in an intermediate critical regime, where 
gauge fluctuations are suppressed, leaving a quenched vector potential. In establishing the limiting magnetic length 
scale we consider the isotropic case, to simplify the notation. In analogy to uniformly rotating supcrfluid 4 He it does 
not depend on temperature [6]. For this reason we set T = T c so that in zero field £ T (T c , H = 0) = oo and for small 
magnetic fields the scaling variable , while for low magnetic fields the scaling variable z = z = £ T _ff 2 /$ tends to 
infinity. In this limit G~ (z) adopts the form [18,34] 

G 3 -(*->oo) = ^ 3 / 2 . (21) 

where C~ is a universal constant. Substitution into Eq.(20) yields for the singular part of the free energy density the 
expression 

= W G - (z) ^ ^Q-C-k B T c ( |-V /2 = Qs ^ = Qz — 2 , (22) 



where 



L,i=.Ja,.U<=[^) ll \ «=(m V \ (23) 
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and a is a universal constant. This relationship between the limiting magnetic length scale and the transverse 
correlation length is also obtained from the scaling form 

(£±(T,H))- 2 = ^ S (z) = ^S-(z), S-(z = 0) = l. (24) 

The divergence of £ T (T c , H) for H — > requires, 

S~ {z -> oo) = az, (25) 

so that 

(f(T ei H)) a = ¥L = lL' H = L%, (26) 

because a/a = 1 follows from Eq.(22), rewritten in the form f s / (ksT c Q^ ) — £ T (T C ,H)~ 3 = L H 3 . Noting that the 
limiting magnetic length scale Lh does not depend on temperature it applies above and below the zero field transition 
temperature. In particular, below T c the transverse correlation length £ T (T,H) is bounded by Lh, 

e(T,H)<L H . (27) 

This implies that in the (H, T)-phase diagram there is no line of continuous phase transitions requiring £ T (T, H) 
to diverge. Furthermore, in the presence of a weak magnetic field the thermodynamic property, O, of a type II 
superconductor takes below T c the scaling form 

Q (t, H) = Q(t,H = 0) Gq (z) ,..z= = . (28) 

It establishes the connection to the scaling form (9) for a system in a confined geometry of length L. 

In practice there are also inhomogeneities of length scale L. As aforementioned, two limiting regimes characterized 
by Lh < L, Lh > L (Eq.(7)) can be distinguished. For Lh < L the magnetic field induced finite size effect sets 
the limiting length scale, because £ T (T,H) is bounded by Lh, while for Lh > L, £ T (T,H) is limited by L, the 
length scale of the inhomogeneities. Since Lh = \f$o/ (aH) (Eq.(23)), both limits are experimentally accessible. 
Lh < L requires sufficiently large and Lh > L small magnetic fields. Moreover, the occurrence of a magnetic field 
induced finite size effect requires that the magnetic field strength exceeds the lower bound H > <J> / (a^ 2 ) (Eq.(8)). 
On the contrary its absence up to the field H provides a lower bound for the length scale of the inhomogeneities 
L > y^o/ (aH) (Eq.(23)). 

Supposing that Lh < L is satisfied, the singularity in the specific heat is replaced by a broad peak reaching its 
maximum value at T p < T c which decreases with increasing field strength. From Eq.(12) we obtain for the shift and 
with Eq.(14) for the height of the peak in the specific heat coefficient the expressions 



l^p I — 



. 



a 



L H \ $0 



H 1 l 2v , (29) 



^ - A ~ M- + B~ _ <S! * ("-EY" 2 " + S -, (30) 



and 



respectively. We note that for a < (Eq.(16)) the peak height decreases monotonically even for small fields. Since for 

fixed magnetic field strength H, \t p \ oc (S.q) 1 ^ , while c(T p ) /T p oc (£q) 3 = V^ 1 , it becomes clear that the critical 
regime considered here can be attained in superconductors with sufficiently small correlation length amplitude only. 
Here|£ p | can be tuned into the experimentally accessible range by increasing the magnetic field, while c(T p ) /T p is 
controlled by the small inverse correlation volume. 

Since the transverse correlation lengths which are parallel to the applied field are limited by the length scale Lh, 
Eq.(29) is readily extended to the anisotropic case for magnetic fields applied parallel to the a, b or c-axis: 
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M= -p 5 - ,I* P I= pH = hrH ( 31 ) 



where 



Lh « aH a > Lh > aH h > L »< ^ (32) 

The specific heat coefficient at T p adopts then the form given by Eq.(5) with the correlation volume appropriate for 
the anisotropic case. This behavior applies when the condition Lh < L (Eq.(7)), extended to the anisotropic case, 

Ljj a < LbL Cl L Hb < L a L Cl L Hc < L a Lb 7 (33) 

is satisfied. Li, i = (a, b, c) denote the length scales of the sample inhomogeneities. Accordingly, the occurrence of 
a magnetic field induced finite size effect requires that the magnetic field strength exceeds the lower bounds 

H a >^-,H b >^-,H c >^-. (34) 
aLbL c aL a L c aL a Lb 

Otherwise the finite size effect is due to inhomogeneities. When the magnetic field induced finite size effect dominates, 
the scaling form (19) of the specific heat coefficient adopts then with Eq.(31) or Eq.(28) the form, 

^-S-W-lH =,(»). , = « M^ ■ (35) 



/ \ aH c ((1, , 

Accordingly, given estimates for A~/a, B~ and determined from experimental zero field data or the magnetic 
field dependence of t p (Eq.(31)), their consistency can be checked in terms of this scaling relation. Consistency is 
achieved whenever the data for the specific heat coefficient in a magnetic field, plotted according to Eq.(35), falls on 
a single curve, which is the scaling function g(y). At T p , where y = y p = —1, we have g(—l) = 1 and Eq.(18) is 
recovered, while for ±y — > it diverges as g (y) oc \y\ a for a < 0. 

At the first-order melting transition there is a discontinuity in the internal energy associated with coexistence of 
the Abrikosov lattice and the vortex liquid. This in turn gives rise to a 5-function peak in the specific heat [37,38] 
and to a singularity in the scaling function G$ (z) of the free energy density (Eq.(20)) at some universal value z m of 
its argument. According to Eq.(20) the melting line T m (Hi) is then fixed by 

On the other hand there is the T p (Hi), where £f (Hi,T) attains the limiting magnetic length scale L^ i . Using 
Eqs.(31) and (32) we find 

MHOI-pg-L) , (37) 

and the ratio 

\t m 

turns out to be universal. Because a limiting length scale, like Lh oc i? -1 / 2 , also leads to a rounding of first-order 
transitions [28], the vortex melting transition is expected to broaden with increasing magnetic field. 
We we are now prepared to analyze experimental data in terms of the finite size scaling theory. 
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III. COMPARISON WITH EXPERIMENT 



We have seen that measurements of thermodynamic properties in an applied magnetic field open a rather direct 
route to trace thermal fluctuations, to estimate various critical properties, to derive an upper bounds for the length 
scales associated with inhomogeneities, to assess the relevance of critical fluctuations and to probe the interplay of 
the finite size effects stemming from sample inhomogeneities and the applied magnetic field. In particular, the rele- 
vance of thermal fluctuations is established whenever there is a remnant of the zero field singularity at T p , shifting 
to a lower value with increasing magnetic field strength. In this case condition (31) is satisfied and the field depen- 
dence of the relative shift t p should evolve according to Eq.(37). However, due to unavoidable inhomogeneities the 
available experimental data do not extend to the asymptotic critical regime, as required to deduce estimates for the 
critical exponents. Since type II superconductors fall in the experimentally accessible critical regime into the 3D-XY 
universality class, we take the critical exponents and critical amplitude combinations, listed in Eq.(16), for granted. 



1460 




76 7a 80 82 84 86 36 90 92 94 96 
T(K) 

FIG. 1. Total specific heat of an untwined YBa2Cu307_5 single crystal with T c — 91.87 K for H c — H\\c. The numbers 
on top of the spike indicate the strength of the applied magnetic field and the location of the melting transition. Taken from 
Schilling et al. [12] 



In Fig.l we displayed the data of Schilling et al. [12] for the temperature dependence of the heat coefficient of a 
YBa2Cu307_,5 single crystal with T c — 91.87 K at various magnetic fields applied parallel to the c-axis (H c ). As a 
remnant of the zero field singularity, there is for fixed field strength a broad peak adopting its maximum at T p which 
is located below T c . As T v approaches T c , the peak becomes sharper with decreasing H and evolves smoothly to the 
zero-field singularity, smeared by the inhomogeneity induced finite size effect. The spike at T m < T p is due to the 
melting transition [12]. Since T p increases systematically with increasing field, condition (33) is satisfied and there 
is a magnetic field induced finite size effect. Accordingly Eq.(37) for the field dependence of the relative shift should 
apply as long as 3D critical fluctuations dominate. In Fig. 2 we displayed t p versus H for the data shown in Fig.l. 
The solid line is Eq.(37) with 

\t p \ = 0.012 H^ 2 " , (39) 

H in T and v given by Eq.(16). In the low field regime, where this asymptotic form applies it describes the data 
remarkably well, while for \t p \ > 0.04 deviations pointing to corrections to scaling become manifest. To fix the 

universal constant a entering Eq.(23) we use for £q c the estimate = \J £ao€bo = V12.4 • 14.7 = 13. 5A, derived from 
magnetic torque measurements [18], Eqs.(37) and (39) then yield for the universal constant a entering the limiting 
length scales the estimate 

a w 3.12. (40) 

In 4 He this value leads with Lq = \Jh] (2am/$l) and — (27r) _ s _1 to w 5 10~ 2 cm, which compares reasonably 
well with 5 10 _2 cm, the estimate derived by Haussmann [6]. Since at the lowest attained magnetic field, H c = IT, 
(see Fig.l) a shift from T c to T p occurs, Eq.(37) yields for the inhomogeneities in the ab-plane the lower bound 
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L ao > 257 A, which is close to the estimates derived from the finite size effect in the zero field specific heat coefficient 
of untwined YBa 2 Cu 3 07_5 single crystals [18,20]. 




H.(T) H e (T) 



FIG. 2. (a) |£ P 1 versus H c for the data shown in Fig.l. The solid line is Eq.(39) with v listed in Eq.(16). (b) Melting line 
Tm (H c ) in terms of |i m |versus H c for the data shown in Fig.l. The solid line is Eq.(41) 

The existence of the first order melting transition of the vortex lattice, seen in Fig.l in terms of the spikes with 
numbers on the top, implies that the scaling function G%(z) in the free energy density (Eq.(20)) exhibits at some 
universal value z m of its argument a smeared singularity. For a magnetic field applied along the c-axis the melting 
line is then given by Eq.(37). In Fig. 2b we compare this prediction with the experimental melting line, derived from 
the data shown in Fig.l, in terms of the solid line, which is 

\t m \ = 0.029 H]' 2 ". (41) 

Together with = 13. 5A, the value used to estimate a (Eq. (??)), we obtain for the universal value of z m the estimate 

z m « 0.1. (42) 

Since a limiting length scale, like Lh oc H^ 1 ^ 2 , leads also to a rounding of first-order transitions [28], the vortex melting 
transition is expected to become broader with increasing magnetic field. This behavior appears to be confirmed by 
the extended study of the vortex melting transition in YBa2Cu307_a single crystals of Roulin et al. [13]. 
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To explore variations between data and single crystals of different provenance we show in Fig. 3 t p versus H ab and 
H c for YBa2Cu3C>6.9 single crystal with T c = 92.56K derived from the data of Roulin et al. [13]. The solid lines are 

\tp (H ab )\ = 0.0033 Hl' b 2v , \t p (H c )\ = 0.0132 U X J 2 \ (43) 
corresponding to Eq.(37) with H in T and v given by Eq.(16). They yield with a = 3.12 (Eq.(40)) the estimates 

£<L = 5.7A, & = 14.4A, V- = (C T Qb ) 2 & = 469A 3 (44) 

Since at the lowest attained fields H c = 0.25T and H ab = 2T a shift from T c to T p is still present, we derive from 
Eq. (34) for the length scales of the sample inhomogeneities the lower bounds 

L ab > 515A, ^L ab L c > 182A. (45) 

The reported specific heat data also revealed melting lines consistent with the scaling form (37), namely 

K (H ab )\ = 0.008 U X J 2 \ \t m (H c )\ = 0.032 R\l 2 \ (46) 

with v = 2/3, which is close to the 3D-XY value (Eq.(16)). Together with the expressions for the t p (Hi) lines (Eq.(43)) 
we obtain for the universal ratio \t p (-Hi)l / \t m = (ozm) 1 ' 2 " (Eq.38)) the estimates, \t p [H ab )\ / \t m (H ab )\ = 0.41 
and \t p (H c )\ / \t m (H c )\ — 0.42, which agree with (azm) 1 ' 2 " = (3.12 0.1) 3 ^ 4 , derived from the independent estimates 
for a (Eq.(40) and z m (Eq.(42)). 




\ ' ' • 1 • ' • 1 0.104-1 — , — i — , — i — , — i — , — i — , — i — , — i — , — 

-0.02 -0.01 0.00 0.01 0.02 _o.-|2 -0.10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 

t/H c 1 ' 2 " t 



FIG. 4. (a): (c (T, H c ) /T - B~) / (A" /a) versus t/Hl /2u for YBa 2 Cu 3 06.9 with T c = 92.56 K derived from the data 
of Roulin et al. [13], with a and v given by Eq.(16), (c(H,t) /T - B~) in J/K 2 gat and H in T. The values for A~ /a and B' 
are listed in Eq.(48). The solid line is Eq.(47); (b): Specific heat coefficient c/T versus t of YBa 2 Cu3C>6.9 with T c = 92.56 K 
at various magnetic fields derived from the data of Roulin et al. [13]. The solid lines are Eq.(13) with a and A + /A~ given by 
Eq.(16) and the parameters listed in Eq.(48). The deviations of the zero field data from these lines reveal around t — the 
finite size effect due to inhomogeneities. 



As the data of Roulin et al. [13] are rather dense and extend close to zero field criticality, a more detailed finite 
size scaling analysis appears to be feasible. Using Eq.(35), it is seen from Fig. 4a that data collapse is obtained over 
a rather wide range of the scaling variable. Indeed, finite size scaling requires the data to collapse on the universal 

$o/ (aH c (£oc) ) ) • Furthermore, scaling also requires that at T p: where y = y p = — 1, 

and for ±y — > 0, g (—1) = 1 and g (y) cx |y|~' Q ' for a < should hold, respectively. Using £j c = 14. 4A (Eq.(44)), it is 

1 /2 

seen from Fig. 4a that at y p , where \t\ /H c = —0.013, the scaling function is close to 1, as required. In addition, the 
solid line which is 

C ^ /T - B --0M7(^LY, (47) 



(A-\t\- a /a) 



hV 2 
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also confirms the divergence of the scaling function at ±y — 0. To qualify the remarkable collapse of the data we note 
that the estimates for A~ = A~ /a and B~ have been derived from the zero field data displayed in Fig. 4b. The solid 
line for t < is Eq.(13) with 

^—=A~= -0.073(J/K 2 gat), B~ = 0.179(J/K 2 gat), B+ = 0.181(J/K 2 gat), T c = 92.56K (48) 

and the critical exponent a listed in Eq.(16), to indicate the inhomogeneity induced deviations from the leading zero 
field critical behavior of perfect YBa2Cu306.g. Since the data collapse onto the universal finite scaling function 
has been achieved without any arbitrary fitting parameter, the consistency with the magnetic field induced finite 
size scenario is well confirmed. Furthermore, Fig. 4b shows how the broad anomaly in the specific heat coefficient 
sharpens and that the maximum height at T p (H c ) increases, evolving smoothly to the zero field peak, rounded by 
inhomogeneities. 

To estimate the correlation volume from the specific heat coefficient in terms of the universal relation Eq.(15) we 
note that A~ = T c a A~, A-(cm- 3 )= (l0 4 /k B /V gat ) A- (mJ/K/cm 3 ) and V gat = 8cm 3 give A~ = 7.96 lO^A" 3 , 
yielding with R~ — 0.815 (Eq.(16)) the correlation volume V~ w 680A 3 , which is reasonably close to the value derived 
from the magnetic field induced finite size effect (Eq.(44)). As in zero field the correlation volume V~ \tpY Zv cannot 
grow beyond the volume Vi of the inhomogeneities, this limiting volume scale is obtained from 

V i = V-\t p \- 3v . (49) 

and t p evaluated in zero field. With V~ = 469A 3 (Eq.(44)) and \t p \ = 0.0025 taken from Fig.4b, we obtain the 
estimates Vi = 7.5 10 7 A 3 and V^ 3 — 422A which is consistent with the aforementioned lower bounds listed in 
Eq.(45). 

1 /2 

Noting that the critical amplitudes, the anisotropy 7 = (£j fc0 /£S)) an d the length scale of the inhomogeneities 
depend on the dopant concentration [19] it is instructive to analyze the data of Junod et al. [17] for undcrdoped 
YBa 2 Cu 3 06.6 with T c = 64. 2K and 7 w 20 [29]. In Fig. 5 we displayed the data for \t p \ versus H c together with 

|t p (ff c )| = 0.045 jry 2 ", (50) 

corresponding to Eq.(37) with H in T and v given by Eq.(16). It yields with a = 3.12 (Eq.(40)) the estimate 

& = 32.6A, (51) 

compared to (J c = 14 .4 A (Eq.(44)) obtained close to optimum doping. Thus, the critical amplitude of the correlation 
length £q c increases in the underdoped regime, reflecting the approach to the 2D quantum superconductor to insulator 
transition in the underdoped limit, where 7 tends to infinity [18,19]. Typically, the leading scaling behavior becomes 
valid when |i| becomes small and H is small. However, one needs to know: How small is small enough? As there are 
corrections to the leading scaling behavior we know, when \t\ is no longer small, the temperature dependence of the 
transverse correlation length scales as [15,16,18] 

£(t)=&\tr(l + aer\tr + ..) (52) 

where ojv is the correction to scaling exponent, adopting in the 3D-XY universality class the value [16] 

ujv = 0.53. (53) 

The dimensionless correction amplitude g^t will in general become larger near a crossover. As YBa 2 Cu307_,5 under- 
goes with reduced oxygen concentration a 3D-2D crossover, a^r is expected to increase, while the critical regime where 
the leading term dominates shrinks. To substantiate this expectation and to illustrate the effect of the corrections to 
scaling, we set £j (T, H) = Ln a , to obtain with Eqs.(31) and (32) for H c (t p ) the expression 

He {tp) = -^-2 K\ 2u (1 + \t p ry 2 . (54) 

A fit to the experimental data shown in Fig. 5b yields 

H c = 33.97 \t p \ 2v (1 - 1.61 \t p rV 2 , (55) 
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with v and uiv given by Eqs.(16) and (53), respectively. This gives £q c — 40.86A which is considerably larger than 
£q c — 32. 6A, the value derived without the correction to scaling. On the contrary, a fit to the data displayed in 
Fig. 3b yields H c — 319.66 |i p | 2iy (1 — 0.014 |£ p | wiy ) 2 , revealing that close to optimum doping the correction is up to 
\t p \ = 0.04 negligibly small, while in the underdoped sample it is significant down to \t p \ = 0.02. This confirms the 
expectation whereupon the 3D-2D crossover in the underdoped regime enhances the correction to scaling and reduces 
the critical regime where the leading scaling term dominates. From the data of Junod et al. [17] we also deduce that 
at the lowest attained field, H c = 0.2T, a shift from T c to T p occurs. This yields with Eq.(34) the lower bound 

yfLM > 576A (56) 
compared to 515 A, the value close to optimum doping (Eq.(45)). 




FIG. 5. (a) \tp\ versus H c for YBa2Cu30e.6 with T c = 64. 2K derived from the data of Junod et al. [17]. The solid line is 
Eq.(47) with v listed in Eq.(16). (b) H c versus |t p | of the same data. The solid curve is Eq.(55) which takes the correction to 
the leading scaling term, H c = 33.97 |t P | 2l/ , displayed by the dashed line, into account. 

Next we turn to Nb77Zr23 [5], 2H-NbSe2 [4] and MgB2 [3], type II superconductors supposed to have comparative 
large correlation volumes. In Fig. 6 we displayed |t p | versus H derived from the respective experimental data. In 
contrast to the corresponding plot for nearly optimally YBa2Cu307_^ (Fig. 3), the data points to a linear relationship. 
Consequently, the critical regime, where \t p \ oc iZ 1 / 2 " with v » 2/3 holds (Eq.(37)), is not attained. Nevertheless, 
there is clear evidence for a magnetic field induced finite size effect, because T p shifts monotonically to a lower value 
with increasing magnetic field. Since the data points to an effective critical exponent v « 1/2 which applies over an 
unexpectedly extended range, we use Eq.(37) with v = 1/2 to derive estimates for the amplitude of the respective 
transverse correlation lengths and the correlation volume in terms of 

l*,l=^ = -f^ = ^|^, (57) 

with a = 3.12 (Eq.(40). The respective straight lines in Fig. 6 are this relation with the parameters listed in Table 
I. For comparison we included the corresponding parameters for YBa2Cu30g.9 and YBa2Cu30g.6 where the critical 
regime is attained. It is evident that Nb77Zr23, 2H-NbSe and MgB2 are type II superconductors with comparatively 
large correlation lengths. Compared to YBa2Cu306.g and YBa2Cu306.6 the correlation volume is 3 orders of mag- 
nitude larger, which renders the amplitude of the specific heat singularity very weak (see Eq.(13). Nevertheless, 
the unambiguous evidence for the magnetic field induced finite size effect reveals that fluctuations, even though not 
critical, are at work. For this reason there is no critical line T C 2 (H) of continuous phase transitions, but a line T p (H) 
where the specific heat peak, broadening and decreasing with increasing field, adopts its maximum value and the 
correlation length attains the limiting magnetic length scale Lh- Because the fluctuations are also subject to the 
finite size effect arising from inhomogeneities with length scale L, the magnetic finite size effect is observable as long 
as L > Lh ■ The resulting lower bounds for the length scale of inhomogeneities affecting the thermodynamic properties 
are also included in Table I. Noting that these bounds stem from studies where no attempt was made to explore the 
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low field behavior, it is conceivable that the actual length scale of the inhomogeneities is much larger. To our best 
knowledge, the only absolute reference stems from the finite size scaling analysis of the zero field specific heat data of 
nearly optimally doped YBa 2 Cu307_,5, where L was found to range from 290 to 419A [20,18]. Interestingly enough, 
the largest bound found here applies to the cubic superconducting alloy Nb77Zr23. In any case, the lower bounds 
ranging from 182 to 814A raise serious doubts on the relevance of the nanoscale spatial variations in the electronic 
characteristics observed in underdoped Bi-2212 with scanning tunnelling microscopy [21-24]. 




H(T) 

FIG. 6. \t p \ versus H for Nb77Zr23 [5], 2H-NbSe2 [4] and MgB2 [3] derived from the respective references. The straight lines 
are Eq.(57) with the parameters listed in Table I. 





r c (K) 


b, bab 


b c 




a (A) 


vr(A J ) 


L, Lafc(A) 


y/L ab L c {A) 


Nb77Zr 2 3 


10.79 


0.1 




55.1 




1.7 10 b 


> 814 




2H-NbSe 2 


7.1 


0.068 


0.22 


42.91 


93.87 


1.7 10 5 


> 258 


> 258 


MgB 2 


35 


0.09 


0.28 


51.73 


110.24 


2.9 10 b 


> 364 


> 258 


YBa 2 Cu 3 6 .95 


92.56 


0.0033 


0.0132 


5.7 


14.4 


4.7 10 2 


> 515 


> 182 


YBa 2 Cu 3 6 . 6 


64.2 




0.045 




32.6 




> 576 





Table I: Summary of the estimates for the transverse correlation length amplitudes, correlation volume V c — 
(£oafc) £oc' l° wer bounds for the length scales L ab and \fT7^L~ c of inhomogeneities, derived from the magnetic field 
induced finite size effect. For Nb77Zr 2 3, 2H-NbSe 2 and MgB 2 we used Eq.(57), describing the intermediate critical 
behavior displayed in Fig. 6. Using the fact that at the lowest attained magnetic fields, H = 0.1T (Nb77Zr 23 ) [5], 
Hat = H c = IT (2H-NbSe 2 ) [4] and H ab = IT, H c = 0.5T (MgB 2 ) [3], a shift from T c to T p was observed, we obtain 
with Eq.(34) the quoted lower bounds for the length scale of the inhomogeneities which affect the thermodynamic 
properties. For YBa 2 Cu307_<5 the estimates stem from the critical behavior shown in Figs. 3 and 5, quoted in Eqs.(43), 
(44), (50) and (51), and the lower bounds are Eqs.(45) and (56). 



The occurrence of fluctuations in type II superconductors with comparatively large correlation lengths also implies 
a rounding of the specific heat peak in zero field, where inhomogeneities set the limiting length scale. Such a rounding 
was found in the specific heat data of polycrystalline Mg n B 2 , obtained with high resolution ac calorimetry [39]. 
Further evidence for the relevance of fluctuations can be deduced from the specific heat coefficient at T p . In Fig. 7 we 
show the experimental data for Nb 77 Zr 23 [5]. For comparison we included Eq.(2) in the form 

5sfflSl..(i-no m 

with the parameters 

c = 18.1mJ/K 2 gat, b = Q.12T~ 1 . (59) 
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Here fluctuations enter via the coefficient b (Eq.(57)). For H > 0.2T the data is consistent with a linear relationship 
and the coefficient b derived from the field dependence of T p (see Table I). The difference between the zero field value 
Ac (if = 0,T C ) /T c c = 21mJ/K 2 gat and the limit Ac (if -> 0,T P ) /T p = c = 18.1mJ/K 2 gat is due to the triangular 
vortex lattice, for which the ratio between these coefficients should be 1.16 [2,5]. 




FIG. 7. Ac (if, T p ) /T p versus H for Nb 77 Zr23. Experimental data taken from Mirmelstein et al. [5]. The solid line is Eq.(2) 
rewritten in the form (58) and the parameters listed in Eq.(59). 



Although we concentrated hitherto on the specific heat, it should be kept in mind that the magnetic field induced 
finite size effect also affects other thermodynamic and the transport properties, including the magnetoconductivity. 
As an example we consider the linear magentoconductivity of a bulk sample. In the isotropic case the fluctuation 
contribution scales in zero field as a oc (£ T ) Z [18], where z is the dynamic critical exponent. Together with the 
evidence for static 3D-XY critical exponents there is mounting evidence for z = 2 [18,30,40,41] deduced from zero field 
data of the fluctuation contribution to the conductivity. In the presence of a magnetic field we obtain from Eq.(28) 
the scaling form 



aoc(e) F{y), y = 



5>n 



(60) 



where F(y) is a universal scaling function of its argument Since in an applied field the correlation length cannot grow 
beyond Lh, there is at T p the residual conductivity 



a (Hp) cx 



(z-l)/2 



where, using Eq.(23), 



rr (rp\ _ *0 , 2 „ _ $0 A T 



(61) 



(62) 



This implies that type II superconductors in a magnetic field do not undergo a phase transition to a state with 
zero resistance. In particular there is no resistive upper critical field if C 2 (T). It is replaced by the line H p (T) 
where the resistivity attains the residual value p(H p ) — 1/er (H p ). A characteristic feature is the positive curvature, 
d 2 H p (T) /dT 2 > 0, for v > 1/2. This anomalous behavior has been found in a variety of cuprate superconductors 
[42,43] and close to T c the data appears to be consistent with 2v w 4/3, the value for 3D-XY critical behavior 
(Eq.(16)). 

To summarize, we have shown that in type II superconductors subject to a magnetic field the correlation lengths 
cannot grow beyond the magnetic length scale Lh ■ Invoking the finite size scaling theory and the scaling form of the 
of the scaling form of the free energy density of an anisotropic type II superconductor close to zero field criticality 
we determined Lh and explored the resulting finite size effect. In contrast to its inhomogencity induced counterpart, 
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arising from inhomogeneities with length scale L, Lh oc H^ 1 / 2 can be varied continuously in terms of the magnetic 
field strength. Thus, as long as L > Lh the magnetic field induced finite size effect is observable and allows to 
assess the importance of fluctuations, to extract critical point properties of the homogeneous system and to derive a 
lower bound for the length scale of inhomogeneities which affect thermodynamic properties. Our analysis of specific 
heat data for under- and optimally doped YBa2Cu307_,5, where the critical regime is attained, revealed remarkable 
agreement with the finite size scaling scenario. On the other hand the existence of the limiting magnetic length scale 
does have even away from the regime where critical fluctuations dominate essential consequences. Since the correlation 
length £ T cannot grow beyond Lh, there is no critical line H C 2 (T) retaining the continuous character of the zero 
field transition. To illustrate this point we considered MgB 2 , 2H-NbSc 2 and Nb77Zr 2 3, type II superconductors with 
comparatively large correlation lengths. Although the analyzed experimental data do not extend to the critical regime, 
we have shown that there is a magnetic field induced finite size effect and with that fluctuations at work. Furthermore, 
lower bounds for the inhomogeneities affecting the thermodynamic properties have been derived. Their length scales 
range from 182 A to 814 A. Since the available data does not extend to low fields, much larger values are conceivable. 
This raises serious doubts on the relevance of the nanoscale spatial variations in the electronic characteristics observed 
with scanning tunnelling microscopy [21-24]. In any case, we have shown that the various predictions we have made 
here are experimentally verifiable. In particular, the analogy with finite size scaling is robust and it should be 
possible to determine in exactly what regions the data collapses onto a one-variable scaling formulation and to what 
extent corrections to scaling can describe any deviations. Furthermore, in spite of the consensus that in type II 
superconductors in an applied magnetic field a phase transition to a state of zero resistance should occur, we have 
shown that the magnetic finite size effect points to the opposite conclusion. This is in agreement with the work 
reported in Refs. [25] and [26]. 

We thank H. Angst, H. Keller and C. Meingast for useful comments and suggestions. 
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